Simultaneous Estimation of Dimension, States and Measurements: 

Rigidity Considerations 
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The Gram matrix of the experimentally prepared states and the performed measurements deter- 
mines the associated density matrices and POVM (projective operator valued measure) elements 
uniquely up to simultaneous rotations of all density matrices and POVM elements. In [Stark, 
arXiv:1209.5737 (2012)] we showed how the state-measurement Gram matrix can be estimated 
efficiently, and we have seen that in the generic case involving full-rank density matrices and 
full-rank POVM elements, the measurement data never suffices to uniquely determine the state- 
measurement Gram matrix. In this paper, we first provide a sufficient criterion to test, which 
additional assumptions (such as the projectiveness of some measurements) suffice to guarantee the 
uniqueness of the state-measurement Gram matrix up to discrete symmetries. Subsequently, in- 
spired by the non-uniqueness of the Gram matrix in situations involving full-rank matrices, we 
examine more closely scenarios involving projective, non-degenerate measurements and their pure 
post-measurement states. We find that these experiments can be identified on the sole basis of the 
collected measurement data. Moreover, the Gram matrix is determined uniquely and can be read 
off directly from the measurement data. 



In a recent paper [l[, we introduced a coherent strat- 
egy to simultaneously estimate Hilbert space dimension, 
states and POVMs. More precisely, we showed how one 
can determine the Gram matrix associated with the in- 
volved density matrices and POVM elements. Up to si- 
multaneous rotations of all states and POVM elements 
in the space of Hermitian matrices, this Gram matrix 
uniquely determines the family of prepared states and 
the performed POVMs uniquely. For situations in which 
knowing the states-measurements Gram matrix is not 
sufficient, we have developed a heuristic algorithm Q— 
taking the Gram matrix as input — to find explicit real- 
izations of the performed states and POVM elements. 

In [l|, we observed that generically, the measurement 
data does not uniquely specify the states-measurement 
Gram matrix — a necessary condition for the uniqueness 
(up to simultaneous conjugation with unitary or anti- 
unitary mappings) of the density matrices and POVM 
elements describing the experiment. The purpose of this 
paper is the exploration of the set of Gram matrices 
which are compatible with the experimentally collected 
measurement data. More precisely, after having intro- 
duced the setting in Section [I] we describe in Section Hfl 
a method to test whether specific assumptions (such as 
the projectiveness of some measurements or the purity of 
some states) are sufficient to guarantee the uniqueness of 
the estimated states-measurement Gram matrix G (up to 
discrete symmetries) . Subsequently, in Section lH"fl we de- 
scribe how the outcome of this test can be predicted with 
probability one before we have even computed the full 
states-measurements Gram matrix. This allows for the 
computation of phase diagrams that are universally ap- 
plicable. In Section [TV] we provide some simple examples. 
In [l[ we addressed problems that might arise from miss- 
ing characterizations of the set of Gram matrices that can 
be generated via density matrices and POVM elements. 
In Section [Vj making use of the uniqueness test, we add 
a few remarks to this discussion. 



The counterexample in [l[ showing that generically, the 
measurement data alone does not uniquely specify the 
Gram matrix, is not applicable to situations in case of 
which the family of prepared states as well as the fam- 
ily of performed measurements contain matrices which 
are rank deficient. This inspired the closer investigation 
of data tables associated with pure states and projective 
measurements in Section IVIl We show that data tables 
which are generated by projective, non-degenerate mea- 
surements and their associated post-measurement states 
can be recognized as being generated in this manner. In 
these cases, the states-measurement Gram matrix can be 
directly read off from the measurement data. No compu- 
tations are necessary and the Gram matrix is uniquely 
determined by the measurement data. 

Other works associated with full quantum modeling 
going beyond state or measurement tomography focused 
on the falsification of incorrect assumptions about the 
performed measurements [||, temporal drifts in the de- 
vices and consequences of inappropriate quantum 
modeli ng [5| -t7i • On the side of dimension estimation, the 
works [8l4lC| investigate linear dimension witnesses to es- 
timate the dimension of the Hilbert space associated with 
the performed experiment. 



I. ESTIMATING EFFECTIVE MODELS VIA 
CONVEX RELAXATION 

In this section, we are introducing the notation and 
sketch results from [l[ about the assumption-free infer- 
ence of quantum models via state-measurement Gram 
matrices. To construct quantum models solely based on 
the measurement data, we are guided by Occam's ra- 
zor, i.e., we are aiming at finding the simplest quan- 
tum model describing the measurement data {f w ,vk)w,v,k- 
Here, f w , v k denotes the frequency of measuring outcome 
'fc' in case we have prepared an unknown state p w and 
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measured an unknown POVM (E v k)k- More precisely, 
among all the quantum models which are compatible 
with the measurement data, we search for those quan- 
tum models whose Hilbert space dimension is minimal. 
We are focusing on the asymptotic regime (i.e., the num- 
ber of independently repeated measurements goes to in- 
finity), so that according to Born's rule 

Expanding all the unknown states and POVM elements 
with respect to an orthonormal basis in the space of Hcr- 
mitian matrices, we can express all states and POVM 
elements in terms of real (the space of Hermitian ma- 
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trices forms a real vector space) vectors in Mr . These 
vectors can be concatenated into the matrix 

P=(pi\---\pw\E 11 \--- 1 E 1K | • • • | E V1 1 • • • | E V k) 

(1) 



and 



G = P P. 



(2) 



G is the Gram matrix associated with p\, Evk- Here, 
W and V denote the total number of prepared states and 
the total number of performed measurements. To keep 
the notation simple, we assume that each measurement 
has K outcomes. This assumption leaves the explana- 
tions unchanged. We observe that the measurement data 
appears as off-diagonal block in G: 



(3) 



where 
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(4) 



V /w,(l,l) ■ ' • fw,(V,K) 



Since rank(G) = rank(P), finding the minimal quantum 
model with dimension d is equivalent to finding a Gram 
matrix G of minimal rank which is compatible (in the 
sense of |5])) with the measurement data. More precisely, 
the Gram matrix G associated with the simplest quantum 
model is a solution of the optimization problem 



argmin rank G 
subject to G e Gqm, 

G satisfies §5§ for given V. 



(5) 



Here, we use Qqm to denote the set of Gram matrices 
that can be generated via density matrices and POVM 
elements. Unfortunately, the rank function is non-convex 
and consequently, the optimization problem (O is non- 
convex. Thus, it suffers from local minima. It has even 
been shown that rank minimization is NP-hard [TTj . To 
approximately solve problem (O, one can replace the 



non-convex objective function by its so called 'convex en- 
velope' which is — roughly speaking — the 'closest' convex 
function. This procedure is called convex relaxation and 
it maps the original, uncomputable optimization prob- 
lem to the 'closest' convex optimization problem. In 
we showed that the trace of positive semidefinite matri- 
ces is (up to an irrelevant constant factor), the convex 
envelope of their rank. Hence, convex relaxation of the 
uncomputable problem ((5|) leads to 



argmin tr G 

subject to G > 0, ||G|| < M QM , 

G satisfies © for given T>, 



(0) 



(the bound Mqm plays a minor role). Up to simultaneous 
Hilbcrt-Schmidt rotations of all states and POVM ele- 
ments, G determines the states and the POVM elements 
uniquely. However, states-measurement configurations 
computed in this way do not necessarily correspond to 
positive semidefinite matrices, even though all pairwisc 
inner products satisfy all constraints dictated by G. To 
find explicit, quantum mechanically valid density matri- 
ces and POVM elements, we have developed a heuristic 
algorithm 0] that takes the Gram matrix G as input 
and searches for positive semidefinite matrices satisfying 
all the G-constraints. After having mapped the original 
rank minimization fljjjj onto its convex relaxation ([6]), we 
have addressed in |l| the question whether or not the 
Gram matrices ([3]) are uniquely determined by the mea- 
surement data T>. A simple counterexample showed that 
this is not the case if all the involved states and POVM 
elements arc not rank-deficient. 



II. FINDING SUFFICIENT CONDITIONS FOR 
UNIQUENESS 

We have seen that data T> generated by full-rank 
states and measurements does not uniquely specify the 
states-POVM-configration. Thus, if we are only inter- 
ested in unique solutions, we need to make additional 
assumptions about the states-measurement-configuration 
to guarantee uniqueness. For instance, one could have 
reasons to assume that some of the measurements 
are projective or that some of the states are pure. 
Which assumptions guarantee uniqueness of the states- 
measurements Gram matrix? Answering this question 
is the goal of this section. Recall that G is the states- 
measurement Gram matrix. Assumptions like the pro- 
jectiveness of measurements can be described by fixing 
some entries of G. For instance, assume that the mea- 
surement 'V is projective. Then POVM elements E v k 
and E v i, k ^ /, are orthogonal and the G-entry corre- 
sponding to the inner product between E v k and E v \ has 
to be set to zero. Define il to be the set of G-indices 
that are assumed to be known. Thus, Go represents our 
knowledge about G. 

As mentioned earlier (cf. Eq. ([I])), we can rewrite the 
density matrices and the POVM elements in terms of 
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vectors that themselves build up a matrix P such that 
G = P T P. When is G uniquely determined by our knowl- 
edge Go? We are discussing this question regarding the 
states and the POVM elements as arbitrary Hermitian 
matrices, thus ignoring for example the positivity con- 
straint of quantum mechanics. These constraints nar- 
row down the set of possible states and measurement 
patterns. Consequently, uniqueness without the consid- 
eration of the quantum mechanical constraints implies 
uniqueness under the consideration of the quantum me- 
chanical constraints. We can thus view the column vec- 
tors of the states-measurement matrix 



P 



[Pl-Pn\ 



as arbitrary vectors in M. R , with 

R = rank(G). 



(7) 



(8) 



In Eq. ([8]), we implicitly assume that n > R because 
G = P T P and therefore, rank(G) = rank(P). How- 
ever, demanding that n > R is clearly necessary, be- 
cause otherwise the columns of P could not even span 
the column vector space R fl ' associated with the Hilbcrt- 
Schmidt space carrying the density matrices and POVM 
elements. Moreover, we assume that P has been drawn 
from a non-singular probability distribution. The fol- 
lowing discussion cannot answer the uniqueness question 
fully but it provides a sufficient condition for the exclu- 
sion of smooth deformations of those states- measurement 
configuration P which are compatible with our knowl- 
edge Go. Hence, we are not able to exclude that the 
set of P-compatible Gram matrices is a set of discrete 
points. The existence of smooth deformations is equiv- 
alent to the existence of an 77-dimensional smooth sub- 
manifold TV C R jRxn (77 > 1) in an open neighborhood of 
P, P G TV, with the defining property that 



sr - 
Qk H 







(9) 



for all Q = [q\, q n ] G TV and (k, I) G £1. The dimension 
of TV given by 

dim(TV) =Rn- rank (v f {kl - h) (P), V/( fc i"i ^(pjj 

(10) 

The matrix in the argument on the LHS is called com- 
pletion matrix [12j |. We are going to abbreviate it with 
Cq(P). The following discussion follows remarks of 
A. Singer and M. Cucuringu in [l2| and aims at adapting 
a discussion conducted in 13j to our setting. Note that TV 
is never zero-dimensional because the left-multiplication 
of P with elements of SO(R) (i.e., simultaneous rotation 
of all the columns in P; relative inner products are un- 
changed) produces new and valid configurations out of 
P. Let TV C TV denote the submanifold generated by 
these simultaneous rotations of the columns of P. Thus, 
the existence of non-trivial smooth deformations becomes 
equivalent to TV 7^ TV. Let (<Pn,Un) and ((p N ,U N ) 
denote the local charts of TV and TV with ip'^-(P) G 



U N C R^W, ^(P) G U N C R dim W C R d ^W 
and tp N (P) = tp~^ (P), when appropriately embedding 

R dim(JV) mto R dim(J\0. Let ^ • Kdim(W) R dim(W) dc _ 

note the orthogonal projection onto span(LT^). Since Un 
and U N are open neighborhoods of (p N (P) there exists 
an e > such that 

n B e (^(p))) = u„ n B e (^(P)). 

The existence of smooth deformations is equivalent to 
the existence of points x G (Un H B E (ipJ J 1 (P)))) \ 
(Ufr n B E (if] v 1 (P))) . Such vectors x exist if and only 
if there are elements in Un H B E (ipJ f 1 (P)) such that their 
orthogonal projection onto kei(7Tj^{UN) is non-vanishing. 
This on the other hand is equivalent to dim([/jv) > 
dim(U N ) which is equivalent to dim(TV) > dim(TV). We 
conclude that the existence of non-trivial smooth defor- 
mations is equivalent to dim (TV) > dim(TV). Eq. iflt)]) 
provides a formula for the dimension of TV. In order to 
check whether or not dim(TV) > dim(TV), we are thus left 
with the task to determine the dimension of TV. 

The dimension of TV is equal to the dimension of its 
tangent space TpN at P, 

TpTV = span{Q(0) | Q{t) = cxp(aft) • f)P}, 

with a(-) C M dim ( s °( fl »,a(0) = arbitrary. Here, / car- 
ries all the generators of SO(R). Rewriting P = USV*, 



S 



\ 



Sr 



r,Rxn 



(11) 











U G O(R), V G 0(n), in terms of its singular value 
decomposition and using that the product of orthogonal 
matrices is an orthogonal matrix, we arrive at 

T P N = span{Vy(0) | W(t) = exp(/3(t) ■ I)SV*}, 

with /?(•) C R dim ( s °( fl )),/3(0) = arbitrary. Taking the 

derivative, we get W(0) = 0(0)- I)P. The matrix /3(0)-J 

is an arbitrary skew-symmetric matrix because /?(0) G 
R dim(so(i?)) can ar bitrary. Denote the linear subspace 
of skew-symmetric (R x P)-matrices by A C M. RxR . We 
conclude that 

dim(TV) = dim ASV* . 

The right-action of S stretches A and embeds it into 
R- Rxn . The action of V* then simply turns the rows of 
the resulting matrices according to V. Consequently, 

dim(TV) = dim{ASV*) = dim(A) = -R(R - 1). 

We conclude that the absence of non-trivial smooth de- 
formations is equivalent to dim(TV) = dim(TV) which is 
itself equivalent to 

rank (Cq(P)) = Rn - ^R(R - 1). (12) 
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III. A PRIORI CHECK OF SMOOTH 
FLEXIBILITY 

Let C := R^( w +v K ) be the space carrying all states 
and measurement configurations and set k := Rn — 
R(R - l)/2 (the RHS of Eq. (p])). Following a proof 
sketch in [l2j], we are going to show next how Eq. (fT2")) 
can be tested before we have even computed any states- 
measurement configuration P that is compatible with the 
constraints Gq. More precisely, we have the following 
procedure in mind (see Algorithm [T]): Fix Q and draw 
Q G C at random with respect to a non-singular proba- 
bility measure on C. We claim that then, with probability 
one, rank (Cq(P)) = k if and only if rank(Cn(Q)) = 
whenever P 6 C is drawn from an arbitrary non-singular 
probability measure. Thus, by testing (fT2"j) on a randomly 
sampled Q, we can test the flexibility of P without even 
knowing P and Gq. The knowledge about ft alone is 
sufficient. 

We start by showing that rank(P) = R (recall our as- 
sumption n > R) on almost all P G C. Let P = USV* 
be the singular value decomposition of P. If we restrict 
ourselves to matrices P with the property that the num- 
ber of non-zero singular values appearing in S is smaller 
than R then we are loosing degrees of freedom. Thus, 
we are constrained to a lower-dimensional submanifold 
embedded into C. Consequently, almost all elements in C 
have full rank. So, both P and Q from before have full 
rank and we know the value of k a priori with probability 
one. 

Next, we are proving the asserted relation between the 
flexibility of P and the flexibility of Q by showing that 
the answer to the question whether or not condition (jT2j) 
is fulfilled is the same for almost all matrices in C. Let 

M :={Q GC|rank(C n (Q)) = /«} 

Proving the claim is equivalent to showing that Af is ei- 
ther equal to C or it is a zero-set in C. Define 

w„(Q):= ( dctA ) 2 

ACCn(Q) 

AeC KXK 

(A C Cq(Q) means that A is a submatrix of Cq(Q)). We 
arrive at the auxiliary claim that 

M = {Q 6 C | u> K (Q) ^ 0} (13) 

"C": Since rank(Gsi(Q)) = k there exists a non-zero mi- 
nor of size k. 

"D": By assumption there exist at least one non-zero 
Co(Q)-minor of size k. Thus, rank(Co(Q)) > n. It re- 
mains to show that rank(Cfj(Q)) > n is impossible. This 
is indeed the case because rank(Gn(Q)) > K would imply 
that 

Rn — k = dim(A) 

< dim(A) = Rn - rank(CVi(<3)) < Rn - k. (14) 



This is the desired contradiction. 

From Eq. (|T3")) we conclude that C\Af is the set of 
roots of a polynomial in the entries of Q G C. The 
polynomial is determined by its values on a finite set 
of points. Thus, whenever C\Af is not a zero-set, uj k (Q) 
agrees with the zero function on a set with accumulation 
point, and therefore, lo k (Q) is equal to the zero func- 
tion. But then, C\Af = C. Consequently, cither C\Af 
is a zero set, or C\Af = C. Thus, we have shown that 
by testing p2[) on a randomly sampled Q, we can test 
the flexibility of P without knowing neither P nor Gq; 
knowing £1 is sufficient. We arrive at Algorithm [T] for 
the a priori test of flexibility. Recall that our analysis 
allowed arbitrary matrices P G M. Rxn and did not take 
into account quantum mechanical constraints such as the 
positive-scmidefinitcness of density matrices and POVM 
elements. Such constraints can constrain the flexibility 
of P. Therefore, Algorithm [T] docs not permit any con- 
clusion about flexibility whenever the condition (fT2"|) is 
not fulfilled. Algorithm Q] provides a sufficient and not a 
necessary criterion for the exclusion of smooth deforma- 
bility. 

Algorithm 1 (Cucuringu, Singer 12]) 

Require: fl, R = rank(G), n = W + VK 
1: Draw Q at random from a non-singular distribution. 
2: Compute the completion matrix Cn(<2). 
3: If rank(Cn(Q)) = Rn~\R(R- 1) then, with probability 
one, almost all P G C are not deformable. 



IV. RIGIDITY PHASE DIAGRAMS 

For each tupel (W, V, K, fl, R), Algorithm[T]detcrmincs 
the property 'locally rigid' or 'locally flexible'. The re- 
sults Algorithm Q] produces when varying W and V, can 
be merged to form a rigidity phase diagram (see Fig- 
ure [T]). In the following, when drawing rigidity diagrams, 
we are each time setting R (equal to d 2 ) to a specific value 
and start to vary the number of states W and the num- 
ber of measurements V. Points in these diagrams that 
are incompatible with the constraint n > R, R = d 2 , 
are automatically assigned to the deformable phase; n 
denotes the number of columns of P from Eq. ([7]). 

We have already observed that V by itself does not 
lead to rigidity. To achieve rigid states-measurement con- 
figurations we have to impose more constraints on G. 
Taking seriously the outstanding role of energy measure- 
ments in experiments, we might promote these measure- 
ments to be projective and non-degenerate. But we did 
not manage to find a rigid phase given this assumption 
(e.g., for dimH = 2, the point (W, V, K) = (100,100,2) 
is still in the non-rigid phase). Each projective, non- 
degenerate measurement provides — via the associated 
post-measurement states — a set of dimH orthonormal 
pure states. But even when postulating the energy post- 
measurement states to be pure and orthogonal, we failed 
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FIG. 1: Rigidity phase diagram in case of dim'H = 2 and pure 
states. 



FIG. 3: Rigidity phase diagram in case of dim'H = 2, 
states and projective, non-degenerate measurements. 
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FIG. 2: Rigidity phase diagram in case of dim"H = 2 and 
projective measurements. 



FIG. 4: Rigidity phase diagram in case of dim'H = 3, 
states and projective, non-degenerate measurements. 



pure 



finding a rigid phase. As a further tightening, we try out 
the assumption that all involved states are pure (now let- 
ting loose of the projective energy measurement). Fig- 
ure Q] shows the rigidity phase diagram in case of two 
dimensional systems. 

We conclude that there exists a rigid phase in case 
the experimentalist knows the Hilbert-Schmidt norms of 
the measured states. Figure [2] shows the rigidity phase 
diagram in case one correctly postulates that all the mea- 
surements are projective and non-degenerate. This con- 
dition is even stronger than the condition that all states 
are pure because the projectiveness of non-degenerate 
measurements implies orthogonality of the POVM ele- 
ments. Thus, in addition to Hilbert-Schmidt norms of the 
POVM elements (diagonal elements in G), we also know 
those off-diagonal elements of G that describe the inner 
product between different elements of the same POVM. 



Figures[3J[Hand[5]display rigidity phase diagrams when 
postulating the combination of the preceding two scenar- 
ios, i.e., pure states and projective, non-degnerate mea- 
surements. 



V. ISSUES IN CONNECTION WITH 

9qm + S+ n B||.|| op <Af QM 

To simultaneously estimate dimension, states and mea- 
surements we first estimate the state-measurement Gram 
matrix via the optimization problem ([6]). Here, we have 
to include the entire knowledge Gq as constraints. Af- 
terwards, we compute the rank of G and check via Al- 
gorithm [T] whether or not G lies within the rigid phase. 
However, we are not able to check whether or not G G 
Gqm (cf. Eq. (JSJ)). If we believe that the data was gener- 



6 



□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□□ 
□□□□□ 
□□□□□ 
□□□□□ 
□□□□□ 



□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 



□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 
□□□□ 



□□□□□□□ 
□□□□□□□ 
□□□□□□□ 
□□□□□□□ 

□□□□□□□ 
□□□□□□□ 
□□□□□□□ 
□□□□□□□ 

□□□□□□□ 
□□□□□□□ 
□□□□□□□ 
□□□□□□□ 

□□□□□no- 



lo 15 20 

measurements 



25 



30 



FIG. 5: Rigidity phase diagram in case of dim'H = 4, 
states and projective, non-degenerate measurements. 
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ated by a quantum model with Hilbert space dimension 
equal to y / rank(G) then G G Gqm whenever G is in the 
rigid phase because there is simply no other constraint- 
compatible Gram matrix with that rank (when ignoring 
hypothetical discrete symmetries). However, as long as 
we are missing a precise characterization of Qqm or the di- 
mension of the underlying Hilbert space, we do not know 
how to guarantee easily that the simplest quantum Gram 
matrix compatible with (|26p does not have higher rank 
then the computed matrix G. This possibility can be ex- 
cluded when we are able to find explicit density matrices 
and POVM elements reproducing G. If we manage to 
find explicit density matrices and POVM elements, then 
we are sure that G corresponds to a quantum model, 
the rigidity considerations apply and we can be sure to 
have found the unique simplest quantum Gram matrix 
explaining the data. A heuristic method (taking G as 
input) to compute explicit density matrices and POVM 
elements on the basis of a Gram matrix will be presented 
in @. 



VI. CERTIFIED UNIQUENESS IN CASE OF 
PROJECTIVE, NON-DEGENERATE 
MEASUREMENTS 



Recall the 2x2 block matrix structure 



G = 




(15) 



of the states- measurements Gram matrix G (cf. Eq.Q). 
Assume that the measurement 't' described by the 
POVM (Etk)f=i i s projective and non-degenerate. Thus, 
it is described by d orthogonal rank-one projectors. The 
post-measurement states (ptw)w=i associated with the 
measurement (Etk)t—i are ^ orthonormal pure states. 
Let T>u denote the submatrix of V associated with the 



post-measurement states (ptw)t,=i an d the projective 
measurement (E t k)i =1 - We conclude that 

P D = l d . (16) 

Eq. (fT6)l forms a necessary condition for (E t k)k=i to be 
projective and non-degenerate. Let Gn denote the (2d x 
2d)-submatrix of G associated with Prj, so that 



Gr 



(G rt )p Vu 
(G m )n 



□ 



> 0. 



(17) 



Next, we are going to show that Eq. (fT5|) forms a neces- 
sary and sufficient condition for (E t k)f = i to be projective 
and non-degenerate, and for (ptw)t>=i to be the associ- 
ated pure post-measurement states. By Cauchy-Schwarz 
inequality, 



INI! \\E t3 \\ 2 2 >\tv(E tjPj )\ = l. 



(18) 



Among all states, the pure states have maximal Hilbcrt- 
Schmidt norm, so that 



lftlli<l- 



Eq. (fT5|) implies that 



\E 



tj\\2 



> 1. 



(19) 



(20) 



On the other hand, we know that (Etk)t = i is a POVM. 
Therefore, V . E t j = I and Ej > 0. It follows that 



d = 



tj 



33=1 (21) 



J' =1 oiH 

Recall that tr(MJV) > for all M,N > 0. Using this 
and Eq. (HO} in Eq. ([2"T]). we conclude that 



\\E t] \\i = i, 

tr(E tj E t] )=0, for all j ^j, 

and therefore, 

((G m ) n ) fj = 5 fj 
for all This in turn implies (cf. Eq. (fT8 
||p«||2>l=HIP«lla = l. 

and therefore, 

((G st ) n ) w = 1 



(22) 



(23) 



(24) 



(25) 



for all j G {1, ...,d}. Consequently, Eq. (fTB)) forms a 
necessary and sufficient condition for {E t k)f =1 to be pro- 
jective and non-degenerate. 
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Assume that all the considered measurements satisfy 
Eq. (fT6|) . Let V and d denote the number of mea- 
surements and the Hilbert space dimension respectively. 
Then, our knowledge of G G IR 2 ™* 2 ™ takes the form 



G 



( 1 * * 




* ' . * 




* * 1 






Id 




* 


\ 


* 



V 



(26) 



* 



Zooming back into the specific measurement 'f ' from be- 
fore using Eq. (fTO|) . we arrive at 



G r 





* 


* 


> 


* 




* 




* 


* 


1 




Id 


Id, 



\ 



(27) 



It follows that 



1 = tr(p^) HMI2 11^112 = 1 



(28) 



However, equality in Cauchy-Schwarz can only happen 
if the involved vectors are linearly dependent. Therefore 
there exists a scalar k, such that p t j = nE t j. Hermiticity 
of ptj and Etj implies that k is real. Taking the norm 
on both sides of the equation p t j = nE t j, we get |k = 
1. Hence, k = 1 because both p t j and E t j are positive 
semidefinite. We conclude that 



Ptj 



E. 



1 j 



(29) 



(30) 



(31) 



(32) 



showing that T> is the simultaneous Gram matrix of the 
states and the POVM elements. We conclude that, once 
Eq. (|16| is fulfilled for every measurement, then the data 
T> uniquely determines the associated Gram matrix via 
Eq. ((311). 



for all t,j. Therefore, 

tv(p t jE si ) = tt(ptjpai) = tv(E t jE si ) 
implying that 

G s t = D = G m . 

Consequently, 

G = 




VII. CONCLUSIONS 



The freedom of choosing the states-measurements 
Gram matrix G can (up to discrete symmetries) be re- 
garded as flexibility of the matrices generating the Gram 
matrix G (under the constraints imposed by our knowl- 
edge of G). Fixing the rank of G and the index set f2 
marking the known G-entries, we showed, using insights 
from rigidity theory, how we can test a Gram matrix G 
for its flexibility before we even know the entries of G 
explicitly. The specification of the index set f2 suffices. 
Determining the rigidity of G for different numbers of 
states, measurements and outcomes, we arrived at rigid- 
ity phase diagrams. To apply this method, we had to 
allow for all G being positive semidefinite instead of G 
being elements of the set C/qm (the set of Gram matrices 
which can be generated quantum mechanically). Since 
all elements in (5qm are automatically positive semidef- 
inite, being in the rigid phase is a sufficient but not a 
necessary condition for rigidity. 

The counterexample in [l[ which showed that V does 
not uniquely determine G applies not to experiments in- 
volving rank-deficient states and POVM elements. This 
lead us to study more closely data tables V associated 
with projective, non-degenerate measurements, and we 
observed that these data tables can be recognized as be- 
ing generated by projective, non-degenerated measure- 
ments. Moreover, G can be directly read off from the 
measurement data. Consequently, the measurement data 
determines G uniquely and no computations are needed 
to determine G. 
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